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It is demonstrated that the oscihations in the width of the momentum distribution 
of atoms moving in a phase-modulated standing hght field, as a function of the 
. modulation amplitude A, are correlated with the variation of the chaotic layer width 
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in energy of an underlying effective pendulum. The maximum effect of dynamical 
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j-H . and minima, respectively, of the chaotic layer width. It is also demonstrated that 
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localization and the nearly perfect delocalization are associated with the maxima 



kinetic energy is conserved as an almost adiabatic invariant at the minima of the 
chaotic layer width, and that the system is accurately described by delta-kicked 
^ ■ rotors at sufficiently large zeros of the Bessel functions Jq (A) and Ji (A). Numerical 
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calculations of kinetic energy and Lyapunov exponents confirm all the theoretical 
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Nonautonomous Hamiltonian systems exhibit a wide range of novel effects with different 



^ . manifestations in the quantum and classical domains [1,2]. A situation of especial interest 
I occurs when the classical dynamics undergoes a transition from stability to chaos as a control 
parameter is varied [3] . In this regard, dynamical localization (DL) is a specifically quantum 
phenomenon taking place in systems where the classical dynamics exhibits chaotic diffusion 
in momentum space, while the asymptotic quantum momentum distribution freezes to a 
steady state although it is initially similar to its classical counterpart. This phenomenon 
is attributed to quantum interferences among the diffusive paths which on average over 
long times are completely destructive. These interferences giving rise to chaotic diffusion 
are properly described in terms of gradual dephasing of the various Floquet states that 
contribute to the initial state. After a characteristic time (the so-called Ehrenfest time 
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te ~ hi. (S/h) /A+, where A+ is the sum of positive Lyapunov exponents (LEs) and S is 
a typical classical action [4]), the momentum distribution covers all significantly populated 
Floquet states as these become completely dephased, leading to the freezing of the diffusive 
growth. This phenomenon has been discussed in model systems such as the kicked rotator [5] 
as well as in the center-of-mass motion of a single ion confined in a Paul trap and interacting 
with a laser field [6] . While no classical or semi-classical explanation for DL has as yet been 
proposed, a formal analogy with the Anderson localization appearing in disordered systems 
has been discussed [7]: DL occurs over the time coordinate with localization in momentum 
space, while Anderson localization refers to spatial coordinates. Additional evidence of the 
DL phenomenon has been provided by the momentum transfer from a modulated standing 
light wave to a sample of ultracold atoms [8-12]. It was found experimentally [9,11] and 
numerically [10,11] that the width of the atomic momentum distribution exhibits oscillations 
as a function of the modulation amplitude, while a theoretical analysis [8] predicted its overall 
monotonous decrease for sufficiently large modulation amplitudes {quantum regime). Also, 
a preliminary qualitative analysis of such oscillatory behavior was previously discussed in 
terms of the zeros of Bessel functions [10,11]. However, no theoretical approach has as yet 
been proposed to explain in detail the waveform of the aforementioned oscillations or the 
amplitude values at which the maximum effect of DL had been found experimentally and 
numerically. In this Letter, such fundamental questions in terms of the strength of homoclinic 
instabilities are explained by studying the dependence on the modulation amplitude of 
the width in energy of the chaotic separatrix layer associated with an underlying effective 
perturbed pendulum. Remarkably, homoclinic instabilities also imprint a clear signature in 
the quantum spectra of a particle confined in a stadium billiard [13]. 

Effective model.— The motion of an atom in a phase- modulated light field, as cre- 
ated for instance by an oscillating mirror, is described [8] by the Hamiltonian H — 
I {2M) — Vo cos [2kx — A sin (cot)], where M, x, p, Vq, k, A, and u denote the mass, position, 
momentum, height of the periodic potential, dimensionless modulation depth, wave num- 
ber, and frequency modulation, respectively. Switching to scaled dimensionless variables 
r = ut.x = 2kx,p = [2k/ (Ma;)]p, H = [4A;^/ {Mu!^)]H, one obtains the dimensionless 
Hamiltonian 

H—— — K cos (x — A sin r) , (1) 
where x, p, r, and k are the dimensionless coordinate, momentum, time, and height of the po- 
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tential, respectively. It is worth noting that the Hamiltonian (1) also describes the dynamics 
of a small dissipationless Josephson junction driven by an external sinusoidal current source 
[14]. To discuss an energy-based characterization of the chaotic dynamics for A > 0, one 
rewrites Eq. (1) as = Hq+Hi where Ho = — Kcosa;, Hi = k [cosx — cos {x — \ sin r)]. 
Now one assumes that H may be regarded as the Hamiltonian of an effective perturbed 
pendulum, with Hi being the Hamiltonian perturbation for A > 0. The change of en- 
ergy of the unperturbed pendulum because of the influence of Hi can be determined from 
dHo/dr — Kp[sm x — sin (x — A sin r)] and hence the change of kinetic energy: d((p^ /2) / dr — 
—Kpsm.{x — Asinr). Since homoclinic dynamics is the origin of the chaotic behavior in 
Hamiltonian systems, it is useful to calculate the change of the averaged energy along the 
separatrix xq^±{t) = ±2 arctan {sinh [y^ (r — tq)]}, Po,± (t) = ±2y^sech [y^ (r — tq)] of 
the unperturbed pendulum Hq: {dHg/dr) = J^^ldHg/dT) dr, where t — Tq gives the 
centre of the soliton [2] . Expanding the term sin {x — X sin r) by using the relationship 
exp (iAsinr) = W (^^''") of the Bessel functions, one straightforwardly finds 

(^) = (^^) ^d{X,K)^f2 K (A, ^) + 2bl (A, k) + 2a„ (A, k) 6„ (A, k)] , 

\ / \ / n=l 

(2) 

where a„ (A, = 167rn^J2n (A) Jq'^ (A) csch (^nn/ \/k\Jo (A)| j / k, 6„ (A, = 

47r|J2„_i(A)/Jo(A)| [l + 2(2n-l)Vl«^o(A)|]sech [(2n - 1) 7r/(2V« | (A)|)] . The 
width function d{X,K) provides an estimate of the width in energy of the chaotic separatrix 
layer. Figure 1(a) clearly shows that the behavior of the width function d (A) correlates 
with those of the classical and quantum mechanical momentum distributions obtained in 
Rcf. [10]. In particular. Fig. 1(b) shows that the maxima of d (A) approximately correspond 
to the amplitude values where a maximum effect of DL appears, while the minima of d (A) 
correspond to the values where classical and quantum distributions approximately coincide 
(i.e., DL disappears). Figure 1 also shows that the function d{X) presents two kinds of 
minima: those exactly occurring at the zeros of Jo (A), where the estimate of the chaotic 
layer width becomes null, and those approximately occurring at the zeros of Ji (A), where 
the chaotic layer exhibits a small but finite width. These two kinds of minima suggest that 
the kinetic energy could be conserved as an almost adiabatic invariant [15] at the zeros {Aq} 
of Jo (A) for a certain set of initial conditions, and at the zeros {A^} of Ji (A) for another 
set of initial conditions. Numerical calculations of the maximal LE (see Figs. 2-4) and the 
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kinetic energy (data not shown) confirm these theoretical predictions. Indeed, system (1) 
presents three LEs (A+,0,A~) which satisfy A+ = — A~. When all three LEs vanish, the 
motion is quasi-periodic and the kinetic energy is conserved in the aforementioned sense, as 
in the examples shown in Figs. 3 and 4. Remarkably, one finds that the function d (A) also 
correlates with A+ and hence with (see Fig. 2). Comparison of Figs. 3 and 4 clearly 
indicates that DL strongly depends on the initial momentum, which should be taken into 
account in experiments [16]. 

Dynamics at the zeros of Jq (A) and Ji (A).— Using well-known properties of J„ (A) and 5t 
(the Dirac comb of period T) [17,2], one straightforwardly obtains that Eq. (1) at {Aq} , {A^} 
reduces to 

l^^+KoST{T")smz = 0{1/K'), (3) 
^ + K^St (r') smx = 0[Jo (At) /XI], (4) 

respectively, where r" = ^/11{t — n/A) , r' = ^/k (t — 7r/2) , Kq = 'k^/k \ J\ (Aq)| /2, K\ = 
ti\/k\Jq (At) I , T = T^\fK, z = X ± tt/2 with the sign + (— ) when Ji (Ag) < (>) 0, and where 
the shift X ^ X ±TT is applied when Jq (^i) < 0. One sees that the description by means of 
the (5-kicked rotors (3), (4) is almost perfect in the quantum regime (i.e., for sufficiently large 
zeros), and that the shift z = x±7r/2 in their spatial coordinates as well as the different 
values of the amplitudes Kq, Ki help to understand why the approximate conservation of the 
kinetic energy does not occur for the same set of initial conditions at both kinds of zeros (cf. 
Figs. 3 and 4). Therefore, Eqs. (3) and (4) can be reduced in the quantum regime to the 
standard maps [2] : p = p + Kq sin z , z = z + Tp; p = p + Ki sin x , x = x + Tp, per period T, 
respectively. In such a case, the energy of a trajectory randomly fiuctuates inside the chaotic 
separatrix layer whose width is AEs^ (Aq, k) « 647rAC~^ | Ji (Aq)|~^ exp — 7r-\/2/ {k | Ji (Aq)!) , 
AEg^i {XI, k) 327rK~^ |Jo (At)p^exp (^—^^/^/K\Jo{Xl)\J, respectively. Note that these es- 
timates are in excellent agreement with Eq. (2). For example, for the experimental value 
K = 0.36 [9,11] one has A^,,o (8.6537,0.36) ~ 0.00138, AE,,i (7.0156,0.36) ~ 0.06574 while 
d (8.6537, 0.36) ~ 0, (7.0156, 0.36) ~ 0.06769. In general, AEs,o (A*, 0.36) , AEs,i (At, 0.36) 
are very small widths at sufficiently large values of Aq, At, which explains the approximate 
conservation of the kinetic energy at such zeros for the majority of the initial conditions. 

Conclusion.— Now an interpretation of the above results in the framework of Floquet 
theory is in order: They indicate that increasing the strength of chaos implies decreasing 
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te, i.e., the time interval in which the Floquet states significantly contributing to the dy- 
namics become completely dephased. Thus, the maxima of the chaotic layer width could be 
associated with corresponding minima of rg. This means that when the width in energy of 
the chaotic layer reaches a maximum as a function of the modulation amplitude, the corre- 
sponding quantum momentum distribution covers a minimal set of significantly populated 
Floquet states, and the effect of DL is maximally enhanced. This picture could well explain 
some previously discussed ways of breaking DL, such as deviation from strict periodicity in 
the driving [18]. Indeed, it has been demonstrated that the suppression of chaos [19] and 
the destruction of DL [18] by a secondary irrational frequency are both very sensitive phe- 
nomena. It is also coherent with the destruction of DL due to dissipation [20]: decreasing 
chaos by increasing dissipation is a universal mechanism [2]. Thus, the present theoretical 
approach to DL suggests that diverse methods for controlling chaos [19] could also be useful 
in controlling DL. Finally, exact experimental realization of the standard 5-kicked rotor by 
using ultracold atoms interacting with a standing wave of laser light can be achieved at the 
aforementioned suitable values of the modulation amplitude. 
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A. Figure Captions 



Figure 1. (a) Comparison between the classical [thick blue (black) line] and quan- 
tum mechanical [red (gray) line] momentum distribution of atoms moving in a phase- 
modulated standing light field for r = 320 (cf. Ref.[10]). Here, Ap = PRMs/i'^f^k) = 
[ip^) ~ (p)^] / (2/iA;) is a normalized RMS momentum width. Also plotted for comparison 
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is the width of the chaotic layer [Eq. (2), thin black line] vs the amplitude A. (b) Differ- 
ence between the classical and quantum mechanical momentum distributions /\pc-Q (green 
line), width of the chaotic layer (black line), and Bessel functions 3 \Ji (A)| (dashed line) and 
3 I Jo (A) I vs the amplitude A. 

Figure 2. Maximal LE corresponding to k = 0.36 and nine initial conditions inside the 
initial well and their average (thick black line), A+, and the rescaled chaotic layer width [Eq. 
(2), red (gray) line] as functions of the amplitude A. The null value of A+ for finite ranges 
of A around the zeros of Ji (A) indicates quasi-periodic motion. 

Figure 3. Grids of 75 x 75 points in the a;(0) — A plane for k = 0.36 and p{0) = 0. 
Blank, yellow, red, and black regions indicate that the maximal LE, A"*", satisfies |A+| < 
10-3, A+ e [10-^ 0.06] , A+ e (0.06,0.12], and A+ > 0.12, respectively. The zeros of Jo (A) 
are also indicated for comparison with the regions where the kinetic energy is approximately 
conserved (blank regions). 

Figure 4. Grids of 75 x 75 points in the a;(0) — A plane for k, = 0.36 and p{0) = — 1. 
Blank, yellow, red, and black regions indicate that the maximal LE, A"*", satisfies lA+l < 
10-3, A+ e [10-3,0.06] ,A+ e (0.06,0.12], and A+ > 0.12, respectively The zeros of Ji (A) 
are also indicated for comparison with the regions where the kinetic energy is approximately 
conserved (blank regions). 
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